Thermodynamics of the anisotropic Heisenberg chain calculated by the density matrix 

renormalization group method 
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The density matrix renormalization group (DMRG) method is applied to the anisotropic Heisenberg 
chain at finite temperatures. The free energy of the system is obtained using the quantum transfer 
matrix which is iteratively enlarged in the imaginary time direction. The magnetic susceptibility 
and the specific heat are calculated down to T = 0.01J and compared with the Bethe ansatz results. 
The agreement including the logarithmic correction in the magnetic susceptibility at the isotropic 
point is fairly good. 



The density matrix renormalization group (DMRG) 
method has brought significant progress to the field of 
numerical calculations. This method provides a way to 
study long chains by iteratively enlarging the system size 
and to obtain the ground state wave function with only 
small systematic errors [jlj . Although the DMRG method 
was initially designed for one-dimensional quantum sys- 
tems it has also been applied to two-dimensional classical 
systems [2], and now even to two-dimensional quantum 
systems [3]. A recent advance in the use of the DMRG 
method is its application to thermodynamics [Q,|5| . In the 
present paper, we apply the DMRG method to the quan- 
tum transfer matrix and obtain accurate results which 
allow us to study low temperature behaviors of thermo- 
dynamic quantities such as the logarithmic correction in 
the magnetic susceptibility of the Heisenberg chain. 

Initially, the DMRG procedure at finite temperature 
was designed to be a similar algorithm to the zero tem- 
perature DMRG except that in the former we take sev- 
eral relevant low-energy states as target states to obtain 
a representation of the Hamiltonian most efficiently at a 
given temperature |l],[| . By adding the density matrices, 
calculated from the low-energy states using a normalized 
Boltzmann weight, one can construct an effective den- 
sity matrix and find important bases by diagonalizing 
it. While this method was successfully applied to finite 
systems, its application to the infinite system is difficult 
because, in general, we need to increase the size of the 
system while decreasing the temperature, which intro- 
duces an enormous truncation error to the calculation. 

To deal with the infinite system in a more system- 
atic way, one can use the transfer matrix method. In 
this method the free energy of the infinite system is 
directly obtained from the maximum eigenvalue of the 
transfer matrix The partition function of the sys- 

tem whose Hamiltonian H is decomposed into two parts 
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where T n is defined by 
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Thus, once we know the maximum eigenvalue A„ 
we can get the partition function 

7 = \ L/2 
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when T n is translationally invariant and the length of the 
system L is large enough. 

The problem now, is how to calculate A max of T, the 
matrix size of which diverges with increasing M . To 
manage this difficulty we use the DMRG method [|: 
We iteratively increase M in T by restricting the num- 
ber of bases used in T and keeping A max optimal. For 
this purpose we define T A (M = 2) = e~^°' iodd and 
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The optimal set of bases at the transformations in eq. (6) 
will be determined as follows. We first diagonalize T = 
7"^ e — /Sofceveng— PohadAfB w j^ n ^he p er i dic boundary con- 
ditions in the j3 direction, and obtain left and right eigen- 
vectors for T, V L and V , which correspond to the max- 
imum eigenvalue of T. We then compose the density 
matrix p A , for the bases of T A e 



as 
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where i and j are indices of the bases for T e 
and e _/3o/lodd T B as shown in Fig. 1. The best way to 
restrict the number of bases is to find a proper transfor- 
mation to new bases i — > i ncw which maximize the norm 
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Y^ij VfiV-j, with a fixed number m of i new . To do this we 
use the left and right eigenvectors of the m largest eigen- 
values j lncw of the pf v to be the transformation matrix 
to the new bases i lew . The error accompanied by this 
truncation is defined by 1 — V™ i 7' new . 
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FIG. 1. (a) Transfer matrix used at the renormalization 
in eq. 6. (M = 2) (b) Transfer matrix used at the renormal- 
ization in eq. 8. (M = 3) 

The most significant difference of the DMRG proce- 
dure for the transfer matrix from the zero-temperature 
DMRG algorithm is that, in the former, we have to use 
the non-symmetric density matrix. This leads to the left 
and right eigenvectors of the density matrix, u* ncw and 
vl ncw , being different. For this reason we have to trans- 
form the transfer matrix using these vectors as 

% A i, ^y"v' ncw T^u^ cw (11) 

ii' 

where t?J ncw and u- ncw satisfy the orthogonality 

E^ neW ^ eW = ^ewW ( 12 ) 

i 

Technically, the non-symmetric nature of the density ma- 
trix is a difficult point of the numerical calculation. In 
fact, if we do not pay attention to the number of states 
kept, the eigenvalue of the density matrix becomes com- 
plex. In order to avoid this we fix the number states in 
the density matrix throughout the calculation. 

In the following we apply the above DMRG method 
to the S=l/2 anisotropic Heisenberg chain to check the 
reliability of the method. The Hamiltonian is given by 

H = JJ2{S*S* +1 + SfSf +1 + A^5f +1 }, (13) 

i 

where A is the anisotropy parameter which changes the 
interaction from XY type to Ising type with increasing 



A: A = corresponds to the XY limit, and A = 1 the 
isotropic point. For this model it is well known that there 
are three different regimes. When A > 1 a finite gap ex- 
ists in the excitation spectrum and both the magnetic 
susceptibility and the specific heat show exponential de- 
cay at low temperatures. At A = 1, the excitation gap 
closes, and the system shows a critical behavior with a 
logarithmic correction in x a t low temperatures ||. In 
the region A < 1, the logarithmic correction in x disap- 
pears and x approaches a finite value in accordance with 
a power law. Thus this model provides a rich variety 
of structures in the thermodynamic quantities and is an 
appropriate test case for the present method. 

In this calculation we use the simple power method to 
obtain the maximum eigenvalue and corresponding eigen- 
vector of the non-symmetric transfer matrix, and use the 
lapack library to diagonalize the non-symmetric density 
matrix. The maximum number of states m in the den- 
sity matrix is 76 and the corresponding truncation error 
is 5 x 10~ 5 at the lowest temperature T = 0.01J. In order 
to carry out the extrapolation for the Trotter number M 
we start the DMRG procedure from several initial tem- 
peratures Pq and iteratively enlarge the transfer matrix 
in the imaginary time direction using the infinite system 
algorithm with decreasing temperature. 

We first show magnetic susceptibility x z obtained by 
the present method. x z is determined from the change 
in the free energy SF — x z h 2 /2 in a small magnetic field 
h = 0.01J along the z-axis. The calculated x z m the XY 
limit is shown in Fig. 2(a). The exact result is also plot- 
ted in the figure |J . The agreement with the exact result 
is excellent even at the lowest temperature T = 0.01J. It 
should be mentioned that such a quadratic behavior in 
X z is difficult to reproduce if we use a simulation based 
on a finite system where % 2 always vanishes exponentially 
or diverges in the low temperature limit. The result at 
A = 0.5 is presented in Fig. 2(b). Similarly to Fig. 2(a), 
X z decreases to a finite value with a power law with de- 
creasing temperature. The obtained power law exponent 
is 2.0 which is consistent with the result of the bosoniza- 
tion ||. The isotropic case is presented in Fig. 2(c). In 
this case the result of the Bethe ansatz is also plotted for 
comparison. The rather rapid decrease in the low tem- 
perature region is due to the logarithmic correction. The 
small discrepancy at T = 0.01J is considered to originate 
from the finite nature of the magnetic field and from trun- 
cation errors in the DMRG. An Ising like case (A = 2.0) 
is presented in Fig. 2(d). In this case x z approaches zero 
with an exponential form consistent with the existence 
of an excitation gap. The estimated value of the gap is 
consistent with the exact value 0.15 J [ jTo| . 
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FIG. 2. Magnetic susceptibility of the S=l/2 Heisenberg 
chain calculated by the DMRG. 

We next calculate the specific heat. The specific heat 
is determined by the second derivative of the free energy; 
C = -Td 2 F/dT 2 . The results for A = 0.0,0.5, 1.0, and 



2.0 are shown in Fig. 3. The exact results in the XY limit 
(A = 0) are also plotted. In the XY-limit the DMRG 
results reproduce the exact values well. For A < 1.0 
the specific heat vanishes linearly as T approaches zero. 
For A = 2.0, on the other hand, it seems to be almost 
zero at finite T consistent with the existence of the spin 
gap. Thus all the results obtained by the present DMRG 
calculation are consistent with the expected behaviors. 
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FIG. 3. Specific heat of the S=l/2 Heisenberg chain cal- 
culated by the DMRG 

To summarize, we have shown that the DMRG method 
at finite temperature is formulated using the transfer ma- 
trix and that this method actually provides reliable re- 
sults even at low temperatures. The magnetic suscepti- 
bility and the specific heat of the Heisenberg chain are 
accurately calculated for both cases in which the exci- 
tation spectrum is gapless and gapfull. The logarithmic 
correction in the magnetic susceptibility is correctly re- 
produced at the isotropic point where the system shows 
the highly critical behavior. Since this method can be 
applied to both fermion systems and frustrated spin sys- 
tems, the DMRG method at finite temperatures will be 
a promising numerical method with which to study the 
low temperature properties of one-dimensional quantum 
systems. 
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